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1. Introduction 

Studies of exclusive semileptonic decays of B, D, and K mesons are used to extract the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements \V u b\, \V c b\, \V CS \, \V c d\, and \V US \ with errors 
competitive with those obtained using inclusive semileptonic decays, leptonic decays, neutrino- 
antineutrino interactions, and z decays [jl|]. The theory inputs needed to fix the CKM matrix ele- 
ments from exclusive semileptonic widths are form factors parameterizing corresponding hadronic 
matrix elements: 

(P2\V^\P l )=f^(q 2 )(p Pl +pp 2 -Ar+f^(q 2 )A^ (1.1) 

where = (m 2 Pi — mpjq^/q 2 , q = p Pl — pp 2 , and V is the appropriate flavor-changing vector 
current. Alternatively we may write [Q] 

where v = ppjmp { and p± = pp 2 — (v ■ pp 2 )v, so that in the rest frame of a heavy meson P\, 



frm-Z&B and 

y/2m Pl v2m Pl p Pl 



Typically, theoretical errors in the form factors limit the accuracy of such extractions of the 
CKM matrix elements. The situation has been acute in the case of D semileptonic decays [Bp. Here 
we describe calculations of D and K semileptonic form factors, which provide access to \V cs u) I and 
\V US \, respectively. 

For D decays we seek not only the CKM matrix elements, but also to validate applying our 
methods to the B decays B — > nlv and B — > Kll. Below we use a subset of the available lattice 
data to check our methods; we compare the shape of a preliminary result for the D — > nlv form 
factor with the shape as measured by CLEO [J3j] . For D decays, unlike B decays, the lattice and 
experimental data overlap throughout most of the q 2 domain, affording a more powerful check. 

Precise determinations of |V MJ | provide stringent tests of first-row unitarity and may furnish 
additional information about the scale of new physics [Q]. Here we describe the main ingredients 
of our strategy to use staggered quarks to obtain ff n (0) and the tests we have performed to verify 
that our approach will yield errors competitive with existing calculations of the form factor. 



2. D -> nlv: Extraction of \V c d\ 

2.1 Ensembles and valence masses 

We have completed generating correlators with Fermilab heavy quarks and asqtad staggered 
light quarks on the 2+1 flavor asqtad staggered MILC ensembles shown in Table [[]. The heavy 
quark is tuned to the charm mass on each ensemble, and the light valence masses include partially 
quenched and full QCD points. In addition to the ensembles shown in Table [T], we are generating 
correlators on a fine ensemble with m/ = 0.15ra v , superfine ensembles with m/ « 0.14m l5 0.1ra 4 , 
and an ultrafine (a ?» 0.045 fm) ensemble with m/ = 0.2m s . However, the analysis presented below 
is restricted to full QCD data from the coarse 0.4m v and 0.2m s ensembles and the fine ensembles 
shown in Table [IJ 
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Table 1: MILC ensembles [^|, |(| [7p for the current round of D — » n{K)lV analyses, together with the valence 
masses used for all ensembles at each lattice spacing. Valence masses after the semicolons are the tuned 
strange mass. Data generation is complete for all ensembles and quark masses shown. 





£3 a (fm) 


ami /am s 


N 3 x N, 


Nconf 


dm valence 


coarse 


0.12 


0.02/0.05 


20 3 x 64 


2052 


0.005, 0.007, 0.01, 






0.01/0.05 


20 3 x 64 


2259 


0.02, 0.03, 0.0415, 






0.007/0.05 


20 3 x 64 


2110 


0.05; 0.0349 






0.005/0.05 


24 3 x 64 


2099 




fine 


0.09 


0.0124/0.031 


28 3 x 96 


1996 


0.0031, 0.0047, 0.0062, 






0.0062/0.031 


28 3 x 96 


1946 


0.0093, 0.0124, 0.031; 






0.0031/0.031 


40 3 x 96 


1015 


0.0261 


superfine 


0.06 


0.0072/0.018 


48 3 x 144 


593 


0.0036, 0.0072, 0.0018, 






0.0036/0.018 


48 3 x 144 


668 


0.0025, 0.0054, 0.0160; 
0.0188 



2.2 Correlators, correlator ratios, and ratio fits 

To extract the matrix elements (n\V^\D) corresponding to /iPj(<Z 2 ), we use ratios of 3-point 
to 2-point correlators designed to cancel oscillations of opposite-parity states in staggered corre- 
lators [§]. To minimize statistical errors and avoid excited-state contamination, we generate the 
3-point correlators at two source-sink separations [^]. The structure of the 3-point correlators is 
shown in Fig. [j]. The 3-point correlators are computed with current insertions at all times between 
the source and sink. For insertion times far from the source and sink, plateaus appear in the ratios. 
These plateaus are proportional to the desired form factors f? n and f® % . 

We are studying combinations of fit ranges, fit functions, source-sink-separations, and mo- 
menta to minimize errors and control excited-state contamination. Fits to ratios yielding fP x on 
the coarse 0.lm s ensemble are shown in Fig. ||; we fit to a constant with an exponential on the 
D-side of the 3-point correlator to account for leading excited-state contributions; the resulting 
curves are consistent with the data. In Fig. || we also plot the resulting plateau terms and bootstrap 



T = 20/21 for p - 0, and T - 16/17 tor p > 



V„ [S] (t) 




D [K] (T) 




p = 000. %/n = 0.76 
p = 100, %fn=Vl 
p= 110.x 2 /n = 0.56 
p= 111. x 2 /» = 0.52 
p = 200, %ln = 0.70 
boot si rap errors 

, , 1 



V4 insertion time t 



Figure 1: Structure of the 3-point functions needed Figure 2: Fits of preliminary results for ffr n to 



to calculate fP^\f^ ]. Light quark propagators are 
generated at t SO urce with local sources [random wall 
sources]. An extended charm [strange] propagator 
is generated at T. 



plateaus and excited-state exponentials. Errors are 
statistical only and were obtained from 500 boot- 
strap ensembles. 
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Z7n = 0.16 



Consistency of lattice and experimental D -> Jt/v form factor shapes 



□ 


0.4m coarse 


□ 


0.2m coarse 




0.4m fine 




0.2m. fine 


o 


0.1m fine 




continuum QCD 




5 ensembles of full QCD lattice data (stal. errors only) 
CLEO-c, PRD (2009). arXiv:0906.2983 



Figure 3: A simultaneous fit to SjPT of all f® n 
data from the indicated ensembles. Errors are sta- 
tistical and were obtained with 500 bootstrap en- 
sembles. The black curve is the continuum result 
at physical quark masses and fiducial energies. 



errors over the entire fit ranges. Consistent with 
separation is optimal for p = 0, while the smaller 



q (GeV ) 

Figure 4: Overlay of the ratio f+ K (q 2 )/ff n {q 2 ) 
from the lattice (red curve and orange error band) 
and CLEO (blue points). The orange error band 
shows statistical lattice errors, and the blue er- 
ror bars, the full experimental errors. At q 2 = 
0.15 GeV 2 , the results agree and the errors vanish 
by definition. 

expectations [^], we find the larger source-sink 
source-sink separation is optimal for |p| > 0. 



2.3 Renormalization and blinding 

We need to renormalize the current. This provides an easy way to do a blind analysis. After 
nonperturbatively renormalizing the quark fields in the current, the remaining lattice artifacts at 
leading order in HQET are perturbatively calculable. A subset of our collaboration is calculating 
this correction, which enters the result as an overall multiplicative factor depending on the ensemble 
and valence masses. By including an offset in this factor, the normalization of the form factors and 
the implied value of \V ct / \ is masked from analysts performing fits, thereby eliminating a potential 
source of bias. 



2.4 Chiral-continuum-energy extra-interpolation 

We fit the form factors f® 7 ^ obtained from the correlator ratios to NLO heavy-meson S^PT [ 1C ] 
and extrapolate the results to the physical quark masses and continuum limit. For the comparison 
below of lattice and experimental results, we also used S^PT to describe the energy dependence 
of the form factor, and we supplement the NLO expressions with NNLO terms analytic in the 
quark masses and lattice spacing. Although we include data from only a subset of the ensembles 



(cf. Sec. 2. 1 ), we perform a simultaneous fit to all data from all ensembles included. We include 



momenta through ^-(1,1,0) and obtain statistical errors by propagating the bootstrap errors from 
the ratio fits. A fit to the data for f® 71 is shown in Fig. ||. The results are stable under variations of 
the prior central values and addition of NNNLO analytic terms. 

2.5 Comparison of lattice and experimental form factor shapes 

To compare our result with experiment, we consider the ratio /f :T (^ 2 )//+ 7r (^ 2 ), where q 2 = 
0.15 GeV 2 is a convenient but otherwise arbitrary reference point. The ratio can 
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be fixed from experiment without the CKM matrix element \V cc i\- Using this ratio to compare the 
shapes of the lattice and experimental results also cancels the blinding factor. The use of this type 



of ratio to compare lattice and experimental results was advocated in []11[|. 

In Fig. we overlay the (preliminary) lattice and (currently final) experimental results for the 
ratio /+ 7t (^ 2 )//+ Jr (^ 2 )- The red curve shows the lattice central value, and the orange error band 
shows the bootstrap errors. The blue data points show the experimental central values, and the blue 
error bars, the statistical and systematic errors from the full covariance matrix [|h. 

2.6 D — > Klv: Summary and next steps 

The shape of our preliminary result for the D — > nlv form factor, obtained from a subset of 
our data, closely matches the shape seen in the CLEO data. This agreement encourages us to apply 
our methods to the calculations of the B semileptonic form factors and related searches for new 
physics. The statistical errors in the lattice form factor at the fiducial value q 2 are about 5%, in 
accord with expectations [g]. We are adding to the analysis partially quenched and full QCD data 
from the remaining two coarse ensembles in Table [T] and the superfine 0.4m s and 0.2m s ensembles. 
Estimates of heavy quark errors, the uncertainty propagated from the D*Dk coupling, and other 
systematics are in progress. 

Finally, we are exploring combining information about the energy-dependence of the form fac- 
tors from the z-expansion with information about the quark mass and lattice spacing dependence 
from S#PT by using S#PT to compute the mass and lattice spacing dependence of the Taylor coef- 
ficients in the z-expansion. This approach furnishes an alternative to S^PT for model-independent, 
simultaneous fits of data at all energies on all ensembles, and is similar to, but distinct from, that 



detailed in [12]. 



3. K — > Ttlv: Exploring methodology to simulate at q 2 = 

One of the most significant systematic errors in traditional lattice analyses of K — > nlv arises 
because correlation functions with periodic boundary conditions do not cover the physical region 
of q 2 , so obtaining ff n (q 2 = 0) and extracting \V US \ from experimental data requires interpolating 
between q 2 max and unphysical values of q 2 . Model dependence is introduced by the choice of inter- 
polating function. We want to eliminate this systematic error by using twisted boundary conditions 



to simulate at q 2 ~ 0. This approach was first suggested in [13] and later exploited with 2+1 flavors 
of domain wall fermions in [ jl4| ] and 2 flavors of twisted mass fermions in Jl5|]. 

The other main component of our analysis is the method developed by the HPQCD Collab- 
oration to study D semileptonic decays [|l6|]. This method is based on the Ward identity relating 
the matrix element of a vector current to that of the corresponding scalar current q^ (k\ Vjj 1 '' \K)Z = 
(m s — m q )(n\S ,at '\K) , with 5 = si, and Z, a lattice renormalization factor for the vector current. 
Using the definition of the form factors in Eq. (1.1) and this identity, one can extract / ( f 7t (^ 2 ) at 
any q 2 by using 



f^(q 2 ) = ^-^{n\S\K)(q 2 ). (3.1) 
m K~ m % 
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The kinematic constraint requires /+(0) = /o(0), so this relation can be used to calculate (0). 
The downside of the method is that it gives no access to the shape of ff n , but that is very well 



known from experiment \m. For more details of this method, see Q16|]. 



3.1 Test run: simulation and fitting details 

The main goal of this test run is a realistic estimate of the statistical errors we could achieve 
and an assessment of how easily we can tune the twisting angles to get values of q 2 close to zero. 
For these tests we used about 500-600 configurations from each of the coarse 0.2m s and 0Am s 
ensembles, and ~ 550 configurations for the fine 0.2m s ensemble. Instead of using the asqtad 
action for the light and strange valence quarks, as in our calculation of the D — > nlv form factor, 



we use the HISQ formulation [|17p, which has better control of discretization effects. 

We generate 3-point correlators as shown in Fig. [I] with a scalar insertion at time t and 2-point 
functions for kaons and pions using both local and random wall sources. The latter sources produce 
results with statistical errors 2-3 times smaller than the former, so in the following discussion we 
consider only the results obtained with random wall sources. We inject momentum in the 2-point 
functions by using twisted boundary conditions to generate (one of) the light propagators. For the 
3-point functions, we inject the momentum p = 0%jL in either the kaon or the pion by choosing 
either 0q = 02 = 0, 1 ^ or O = 1 = 0, 02 ^ 0, respectively (see Fig. [I]). The different external 
momenta and resulting q 2 are shown in Table ||| We have obtained two values of q 2 very close to 
zero by tuning the twisting angle from 2-point correlator fits only. To extract the form factor, we fit 
the 3-point and 2-point correlators together, which gives us slightly different values for q 2 , but still 
close enough to zero to avoid any significant interpolation in q 2 . In fact, the values of ff n {q 2 ~ 0) 
that we obtain from the correlators with external momentum injected in the kaon and the pion agree 
within one sigma. 

We repeat the combined fits using iterative averages of the correlators at different values of t 
and T to suppress oscillations due to opposite-parity states ||8|], including the ground state and first 
oscillating contributions. We use fits to the ground state alone and fits including four exponentials 
to crosscheck the central values and errors. 





|0 2 | 


(m) 2 


1.08 
1.06 








0.0227(3) 


1.04 





0.7295 


0.0011(4) 




0.7295 





0.0153(3) 


& 1 





0.9105 


-0.0109(5) 


0.98 


0.9105 





0.0114(5) 


0.96 


1.2876 





0.0003(3) 


0.94 



Table 2: Simulation values of the twisting angles 
01 and 82, and the corresponding q 2 . Errors are 
statistical. The smallest q 2 available with periodic 
boundary conditions is ~ —0.104, already outside 
the physical region. Lines in bold correspond to 



q 2 ~0. 



• 0.010/0.050 coarse ensemble 
■ 0.0062/0.031 fine ensemble 



(r,q) 



Figure 5: Form factor /q k as a function of momen- 
tum transfer q 2 for the coarse and fine lattice points. 
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3.2 Test results and future plans 

The results for the fQ K (q 2 ) form factor for the different values of q 2 simulated on the 0.2m s 
coarse and fine lattices are shown in Fig. |5| The statistical errors for the two coarse points and 
the fine point with q 2 ~ are about 0.4%. In Fig. |5] one can see that the form factors for a ~ 
0.12, 0.09 fm agree with each other within statistics, suggesting very small discretization effects. 
Similar behavior is observed when comparing results with mi = 0.2m s and mi = 0.4m s on the coarse 
lattices. Such behavior suggests that, after extrapolation to the continuum and the physical sea light 
quark masses, residual effects for those error sources will be negligible. 

Using the full statistics available in these ensembles, around 4 times the number of config- 
urations used here, we expect statistical errors around 0.2 — 0.3%. Since we will eliminate the 
uncertainty due to the q 2 interpolation by simulating at q ~ 0, the only significant remaining error 
besides statistics will be the one associated with the chiral-continuum extrapolation. We plan to do 
this extrapolation using continuum #PT at NNLO and incorporate taste-breaking effects at NLO, 
including data from at least three different lattice spacings and light quark masses down to m s /S. 
Based on the tests described here, we expect our calculation to be competitive with the existing 
state-of-the-art calculations of ff K {0) Q, 0]- 
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